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ФИЗИКО-МАТЕМАТИЧЕСКИЕ НАУКИ 
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В.Л. ЗАКОВОРОТНЫЙ, ФАМ ДИНЬ ТУНГ, НГУЕН СУАН ТЬЕМ 


МОДЕЛИРОВАНИЕ И ИДЕНТИФИКАЦИЯ 
ИНЕРЦИОННЫХ И ДИССИПАТИВНЫХ СВОЙСТВ ПОДСИСТЕМ 
РЕЖУЩЕГО ИНСТРУМЕНТА И ЗАГОТОВКИ ПРИ ТОЧЕНИИ 


Рассматривается вопрос математического моделирования и идентификации параметров динамических под- 
систем инструмента и обрабатываемой заготовки при изучении процесса резания. 
Ключевые слова: математическое моделирование, процесс резания, динамические подсистемы. 


Введение. Моделирование динамических свойств подсистем режущего инструмента и заготовки, 
взаимодействующих в процессе обработки через динамическую связь, раскрывает зависимость 
сил от упругих деформаций инструмента относительно заготовки и технологических режимов [1- 
3]. При изучении устойчивости и методов формирования стационарных многообразий в окрестно- 
сти равновесия достаточно анализировать деформационные смещения в плоскости, нормальной к 
поверхности резания [1, 2]. В этом случае можно воспользоваться уравнением, полученным в 
предположении, что деформации заготовки отсутствуют: 

4?Х 
а? 
где Е ={Е(Х,5,„1,), В (Х ты — вектор-функции динамической характеристики процесса, 





т + их“ +<(Х) = Е(Х,8 5), (1) 


раскрывающие зависимость сил резания от упругих деформационных смещений инструмен- 
та, а также от технологических режимов: величины подачи на оборот 5, и глубины резания 


Е, при заданной скорости; Х = {Л ,Х р — вектор упругих деформационных смещений 


вершины режущего инструмента; (Хх) = т, + (|, (Х) = + (Х)], сх) = [6 .0|, 
5,К =1,2 -— функциональные матрицы инерционных и диссипативных коэффициентов соот- 
ветственно, а также функциональная матрица формирования упругой составляющей сил в 
зависимости от вектора деформационных смещений и технологических режимов. Во всех 
случаях здесь и ниже символ Е — операция транспонирования. 


Моделирование и идентификация параметров инерционных и диссипативных характери- 
стик подсистемы режущего инструмента позволяет наряду с ранее полученными результатами 


создать полные математические модели подсистем без процесса резания: (Хх) = т, + (| 


и 7(Х) = [м (Х )] . Однако построение математических моделей и идентификация параметров в 


равной мере могут быть распространены и на подсистему обрабатываемой заготовки, а также на 
другие технологические операции обработки резанием. Свойства динамической связи процесса 
обработки будут рассмотрены в следующих наших публикациях. 

Динамическая модель подсистемы режущего инструмента. Главное внимание уделим 
инерционным и диссипативным свойствам подсистемы инструмента в вариациях относительно 


Т 
точки равновесия ЛХ” ={Х",Х.} , которая определяется из уравнения (1) для постоянных и 
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неизменных внешних сил, действующих на вершину инструмента. Полагается, что матрица 
с(Х) = [+ (| является известной. 

Тогда для малых значений х={х,,х,}' (Х(@=Х" +х(0 ) уравнение динамики пред- 
ставляется в виде 
ах 
4? 
где т(Х“), #(Х”), с(Х") — линеаризованные матрицы инерционных, диссипативных и упру- 





(Х) “нех = 7, 2) 


гих характеристик; Г = {1 (1), Г, ()}"- внешнее 5-образное возмущение, по реакции на 
которое осуществляется идентификация параметров. 

Исследования колебательных реакций инструмента на внешние б-образные возмущения 
показывают, что они обладают пространственной избирательностью. При этом также можно вы- 
делить две системы координат {Х\,Х,}ЕХ и {1,75} Е У, связанных между собой линейной 
операцией поворота координат на угол а. Тогда уравнение динамики в вариациях относительно 
точки равновесия в пространстве У определяется системой из двух независимых уравнений 
(рис.1) 


р 4 
2 ТА СТ су = Г, (3) 
(у) То 0 (у) пол 0 
где т” = : — матрица инерционных коэффициентов; й“” =| ’ — матрица ко- 
0 > 0 №2 
д _|б 0 т 
эффициентов демпфирования; с” = 0 — матрица жесткости; у=\{у,,у›}? - век- 
С0,2 


тор деформационных смещений; /“ = {/®, {®\' - вектор дополнительных сил, возму- 
щающих движения системы в окрестности равновесия. 





Рис.1. Система координат, в которой отсчитывается упругое деформационное смещение 
вершины режущего инструмента на внешние силы 


Системе (3) соответствует следующее уравнение в пространстве Х’. 


ах ф 
+й—+сх=Х, 4 
ЧЕ? '” 1 (9) 





т 
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: 1 ь 
(то1с05* 9, + то › эт? 0) Ь (то > — тол) т 2 


т ТТЛ 
где т = Г — матрица инерцион- 
: : то т 
Е (то-то 1) т го (то т” а + ть > со3” а) р 
ных коэффициентов; 
Г ь 1 : 
(101С0$* 9 + 2 эт? 0) В (о. 2 — Мол) т 24 в А, 
= = ° | — матрица диссипации; 
Е (10.2 — Мол) го. (пол $” 0+ Йо С03° а) с 
Г . 1 . 
(С01С08°@ + Со,2 т” 00) В (Со,2 — Сол) 24 ел С21 
е=| Г =| * ° | - матрица упругости; 
ь : саг -@ 
Ь (Со2 — бол) 2 (Сол т” 9 + Со 05° а) о. 





х=(х,х,} -— вектор деформационных смещений в окрестности точки равновесия; 


Д={/,№}’ - вектор возмущающих сил. 

Как видно, в пространстве Х координаты состояния являются связанными. Прежде чем 
обсуждать алгоритмы идентификации параметров, проанализируем некоторые важные свойства 
систем (3) и (4). 

1. Матрицы шт, Й и с являются симметричными и положительно определенными. Не- 
трудно показать, что определители от матриц т, Й и с для системы (4) соответственно равны: 


: 1 . 
(то1с05° + то т? 0) Ь (ту.> — т) т 2 
1 = То, [№] = йо лйо,2 ; С] = С01С0,2, 
В (то > — тол) т 2. (тот эт? а + ту. с0$? ©!) 


т.е. их значения в пространствах Х и У одни и те же. Все силы в этом случае обладают потен- 
циальными свойствами, и без процесса резания система имеет равновесие асимптотически устой- 
чивое, что естественно. 


2. Характеристические полиномы Д®(р)=АД(р) систем (3) и (4) равны, так как координа- 


ты рассматриваемых систем связаны линейным преобразованием. Равенство характеристических 
полиномов легко доказать и их прямым вычислением. 
3. Если идентифицированы параметры в пространстве У, т.е. ту. ,пто Йо, 


1о.2Сол»Со.2›@ , ТО элементы матриц т, й, с в пространстве Х определяются из зависимостей 


(4). Справедливо и обратное. Если определены элементы матриц т, Й, с, то можно вычислить 
параметры т, 1, то. , Ио 1> Ио. ,Со1›Со.2.@ ИЗ выражений: 


с соза —5ша| | С1 21 с0$0 зто | _ 
Шо —с0$9, 2 22| |-919а 059 


_ (с1160$* 0 + с> 2 эт? а — 2с1.2 91 0605 9) [0,5(с>2 —с11) 91 20 + с, › 60$ 24] о 

Е [0,5(С>2 — 11) т 20 + с, › 60$ 24] (са эт? 9+ >> с05? + 2с1.2 $1 0, с05 9) | 

5 — (а 05? + й› › эт? а — 2й, › 9 9. с0$ 91) [0,5(й> > — #1) т 20, + й, › с0$29] (6) 
[0,5(й> > — #1) 20+ Й, › с0$ 20] (т т? 9+ й› › 05? + 2, › 91 9, с0$ 41) ' 
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д | (И 0$? @ + т>> т” а — 2т,.› п 060$ 0) [0,5( ть > — та) 20 + ти › с0$ 24] (7) 
т — 
[0,5(т. > — т 1) эт 24 + т, › с0$ 20] (т, | эт? 9+ т> > с05? о + 2т; > эп 00$ 0) : 
1 РВ 1 21 1 2т 
© = —агсо - = —агсю — = —агсе -- й (8) 
2 (с 2 — Ст) 2 (й, › = й т) 7. (т. › № тт) 


В частности, справедливость уравнения (8) может быть критерием справедливости схема- 
тизации, приведенной на рис.1. 
Методика идентификации параметров подсистемы режущего инструмента. Исходной 
информацией для идентификации положены экспериментально определенные матрицы весовых 


функций, рассматривающих преобразование сил Г“ (1) = {Х (#),0}' и (= 10, (0 в 
колебательные ускорения Х (Г) и Х,() и их частотные образы. Здесь Г®(Р) и Г®(® - 
б-образные внешние силовые возмущения. 

Исследование матриц т, й выполняется в два этапа. 

На первом этапе определяются параметры системы в окрестности точки равновесия 
Х* ={Х! > 0,Х; > 0%”. При этом для устранения возможных люфтов в неполностью затяну- 


тых соединениях осуществляется предварительное силовое нагружение системы с последующим 
снятием нагрузки без изменения ее знака. 

Вначале вычислим выражения для изображений по Лапласу для колебательных ускорений 
%, (р) и Х, (р) ‚ обусловленных силой /®(р)={Л(р),0}", 

2 2 2 2 

Р[(т.›р +№)р+с›>)] (р). у —Р[(т ор +ор+с>)Л(р)] 
й 2 = 
д(р) д(р) 


где Д(р)= (тр” +Й.р+ се. (т,,р” И. „РС. )= (тр чар ба — характеристи- 





‚ (9) 





х(р)= 


ческий полином системы (3). Колебательные ускорения Х/(р) и Х,(р), обусловленные 


силой Г“(р) = {0, /,(р)}", 
— Р"[@т,>р° ы р +С12 Л (р) а Р"[(@т, р” + р + си) (р) 
А(р) _ А(р) 


Как показано выше, Д®(р)=Л(р). Поэтому для автоспектров колебательных ускорений, 








%(р)= (10) 


вызванных силой /“(р)={/(р),0}', имеем 


©“ [(с> т т, 26’ )” + (й> >6)* 15 ‚Л (©) 














(0) те 
ня [со — тов) + (и.о Пси - тов" + (п >67] т 
5. = © р —т, 2 ® . + (#,›6) ли — 
7? [(Сол — Ту”) + (16) (Со, — То)" + (й, >60)" ] 
а для автоспектров, обусловленных силой (р) = {0, /,(р)}",- 
5®. (@)= о" __ _ т, 26 й + (зв) бд —: 
—_ [Сол ки ) + (иле) (о о У + (1.262) ] (12) 
5% (®)= © Ко ИН + (1,16) А — 
к. [(Сол — То”) + (16) (С — ть”) + (№26): ] 


При оценивании параметров системы проанализируем выражения (11), (12). Во-первых, 
спектры 5, , (0), 5, ;, (©) и 5, ‚(@) являются измеримыми, причем, спектр 5’, , (6) — широ- 
кополосный, и в рассматриваемом частотном диапазоне его можно считать белым шумом. Широ- 
кополосность спектра обеспечивается подбором бойка в ударном молотке. Во-вторых, структура 
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спектра такова, что он состоит из двух резонансов (©, и ©. ,) и одного антирезонанса <>, , 
(рис.2). Для оценки т, и т,, достаточно определения автоспектров, так как системы являют- 
ся достаточно добротными. Поэтому определение частот (<>, и <>, , ) по автоспектрам не вызы- 


вает сложностей. Известно, что резонансная частота, соответствующая максимальной амплитуде 
в рассматриваемой диссипативной системе, отличается от собственных частот (©, и ©, ,), по- 


лучаемых из условий: 





































































































































2 ы 2 — 
Сил 1 (2. ,) = 0; => (0.2) = 0, (13) 
где о, и Фе — резонансные частоты соответствующей консервативной системы. 
№10 
5 1110* 
5 
Ве.10° 
0 6 Ве- 10° 
Г 
== 
ы. © =1200 с" -5 
-10 -10 
® =2000 с" р. Е 
-15 -15 
Фи =2000 с”! 
-20 -20 
-1 -0.5 0.5 1 1 
а) 
№110 
0.5 —_ 
Ве.10° 
0 
Антирезонанс 
-0.5 
-1 
-1.5 Фи = 2000 с 
-2 
© =1200 с"' 
-2.5 


























-1 -0.5 0 0.5 1 1.5 


в) 





Рис.2. Изменение взаимных спектров колебаний по силам при различных значениях угла В ориентации 


внешнего силового возмущения в пространстве У:а- 5’„, (16), В=0;6- 5. (16), В=л/4; 
в- 5, (10), В=л/4;г- $ у, (10), В=л/2 


Однако это отличие заметно лишь при коэффициентах затухания, больших 0,1. Для опре- 
деления резонансных частот необходимо вычислить экстремумы функций 


[(о1 =. 6”)? + (16) 1 и [(с›-то 26”)? + (>60)? ]. Для первой мультипликативной формы 


(с — то 16°)? Е (116)? имеем 
1 


01° = [(©09.1)? - (40.1)? ]2, 
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5 Сол 2 (о) 
где (О.,) =—^; (Аб) =, 
то 2(то.) 
ИЛИ 
и 
ОР) Е (О.П — 2(Е, )?]2 и 1= 1,2, (14) 


где & - коэффициент затухания рассматриваемого колебательного контура. 


Диапазон значений коэффициента затухания в механических системах металлорежущих 
станков обычно 0,02-0,06. Так как упругие характеристики подсистем известны, то легко опреде- 
ляются обобщенные массы: 


Сол С0,2 
т ПЕ (15) 
(О )" (о) 


Для вычисления значений ти т, не обязательно измерять возмущающие силы, так 


как обобщенные жесткости с.; и с,, известны. Частоты ©, и ©,, определяются на основе 


автоспектров в безразмерном виде. 
При оценивании коэффициентов затухания учтем, что внешние силы измеримы, и их мож- 
но ориентировать по направлениям /^. и Л.,. Поэтому 


5 кк (Эа) 5 к Е (О) 


О (©) 9 


Так как колебательные контуры добротные, то коэффициент затухания можно определить 
также на основе анализа уширения спектральной линии каждой из резонансных частот, т.е. 


Е 












0,2 = 








(16) 


где о. -@0) — ширина резонансной кривой на высоте амплитуд 0,707 от резонансной ам- 


1 
> 


плитуды. 
Тогда коэффициенты демпфирования й,, и й, , равны: 


йо = 25: Ото: = 26, Мо: Со, › 1=1,2. (17) 
В тех случаях, когда резонансные кривые имеют близкие частоты, можно использовать 
выражения, вытекающие из выражений (11) и (12) 
(0 2.2 2 (2) 
ба (0) (с›-т,›®”)” + (№5 56) 5 


А (©) _ 
(0 (2) 
м 5 


№,» 


(©) > (с2-— т, 6”)? + (й,26)* 

ии 22 2" 
(©) (са -тию“)“ +(й 16) 
При этом в качестве опорных удобно взять те частоты, на которых влияние коэффициентов 


демпфирования максимально, т.е. — антирезонансы. 
На втором этапе рассмотрим систему, имеющую произвольную точку равновесия 








А (&) = = . 
(©) (с<2- т, 60°) + (1, 26)* 


Х* ={Х',Х;}". Для смещения точки равновесия необходимо создать внешнее силовое поле. С 


этой целью введем в систему дополнительный упругий элемент (рис.3). Дополнительная жест- 
кость вносит погрешности, прежде всего, при определении собственных частот. Главное требова- 
ние к упругому элементу — линейность упругих свойств. Кроме этого, обеспечим высокую доброт- 
ность изгибных колебаний упругого элемента. Тем самым его диссипацией, по сравнению с дис- 
сипацией подсистемы режущего инструмента, можно пренебречь. Тогда введенная дополнитель- 
ная связь характеризуется исключительно упругими свойствами, определяемыми матрицей 


а к, икра 
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Уточним уравнения динамики с учетом параметров этого элемента. Вместо уравнения (4) 
имеем: 





п ® нео +6. х(0 = ГО, (18) 


ЧЕ? 


), причем, с‘) — матрица жесткости дополнительного упругого элемента. 


где с; =с+с 
В системе координат (Х",Х») матрица с° является диагональной с диагональными 
12 2 


элементами с/у, со,›,. Представим уравнение (18) в системе координат (7, ,7,). 


4 (0) 40) 
т) —+ ро и с?у() = (в й (19) 
[А Ге 
т 0 р 0 
где т“ = о ; 20 = о ; 
0 т 0 йо 
2 2:39 : ^ А 
д _ | бол +Соду 608° @ + су 5 а (Солу — Сор у} 060$ 0) ы о а 
(Солу — Со2у ) 1 9.608 0) С0,2 + Солу т? 9+ Со, с0$? © м а 


Система (19) имеет характеристический полином 

А = (тоар” +Жр+ в тор” +2р+ ера) == а г (20) 
Исходной информацией для оценивания матриц 7, Й являются матрицы с, с”, с; ; угол ори- 
ентации эллипсов жесткости с; виброускорения (х,,Х.) и внешнее б-образное силовое возму- 
щение /(г). Автоспектры колебательных ускорений, определяющие реакции системы на возму- 
щение Г“”(р)={/(р),0}' ‚можно представить в виде: 


в [(с2; тб 50°) + (#59 50) 151 .д (©). 














5%, (©) = 
5: х о) (©) (21) 
50. о: Оо) бо аи 
9, Ао (©) ? 
а автоспектры, обусловленные силой /“(р) = {0, /.(р)}', будут 
5%. (®)= в“ [(с<®; - т. 0”)? +(15.6) 15 5х и 
5, х 44 (©) (22) 
$. (®)= "(с р 5-1? )* Е (11°) 5. (©) 
м А (6) ` 


А 22 2 22 2 2 2 

где ^ (©) = Кат — То ) + (1,16) ее — То 20° ) + (№ 26) 1-2.) Е — То )х 
х (02; —-ть 60°) — И.о. 6)" 1+ (с.)“ - квадрат модуля характеристического полинома 
для р=/ю; с ЕЁ ©], в = орт ® = = О 5,К =1,2 — матрицы жесткости, дис- 


сипации и инерционных коэффициентов в пространстве Х’ соответственно. 
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Эллипс податливости 
дополнительного 
упругого элемента 


Несущая система 
станка 





Рис.3. Схема введения дополнительного упругого элемента, обеспечивающего 
смещение точки равновесия подсистемы режущего инструмента 


Вначале проанализируем функцию Л“ (%) для консервативной системы 
А 22 22 2 2 
д (©) = а — То 16 ) и — То,2 © ) -2(<®;) о — То 16 )х 


2 4 
^ о — То, 26 р в . 


(23) 
Если ориентация дополнительного упругого элемента совпадает с ориентацией одной из осей 
(у,,у>), то в выражении (23) в =0. В этом случае характеристический полином системы сов- 
падает с полиномом без дополнительной связи, т.е. 
А 2\2 2 22 2 
д: ее — т, 60°)” + (#10) а = — ту 26°) + (26) ], 


и определение обобщенных масс и коэффициентов диссипации подсистемы режущего инстру- 
мента осуществляется по отмеченным выше алгоритмам. Однако конструктивное выполнение это- 
го условия вызывает некоторые сложности. 


Проще всего дополнительный упругий элемент ориентировать по направлениям Х| И 
А: 
Вначале рассмотрим алгоритмы оценивания обобщенных масс То И Ту › по уравнению 
консервативной системы. Для этого вычислим экстремальные точки 
(А) =: 2 (^.) (^) =. 
ОА” (6) / дв = 4, (в)[2ту то, >65 — (655 тол +смхто›)]=0, (24) 


22 2 2 
где Л, (©) = (с); —ть 16°) (с ; то?) (с). 


Рассматриваемое выражение имеет четыре экстремальных точки (рис.4): 
(у) (у) 
СУ + СТ , 


(то) + (то)? 





90°) =0; (©®-2))? = 





(у) (7) (у) (7) 2 (7) \2 
11517101 + С2,250,2 — сть № 655 =Ту 2) + 4(с1> ;) 70,170, 


(©, =“ ‚ 025 


2тулто > 








(7) (у) (7) (7) 2 (7) \2 
11.5101 + 62,2. 50,2 + сть — <55 >То2) +4(с15 5) Тодт 


2ту то 





(©, .)*=- (26) 
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Первые две частоты соответствуют максимумам, частоты ©, и <>,, — минимумам, т.е. 
резонансам. При с), =0 обобщенные массы вычисляются из выражений: 
— в) —. — в) —2 
то = 1х (01); тор = 6225 (2) °. (27) 


Вид функций А“ (0) показан на рис.4. 


АА (®) 


А 212 
Е | 





ОР —0 Оо ©.2) Оо [0 


Рис.4. Качественный вид функции 44) (6) 


В зависимости от коэффициентов затухания расположение годографа будет смещаться, 
как показано точечной диаграммой. Однако в рассматриваемом классе систем смещение резо- 
нансных частот будет незначительным, так как добротность рассматриваемых систем достаточно 
высокая. 

Для вычисления обобщенных масс и геометрической интерпретации решения удобно сде- 
лать следующие преобразования. Сложим правые и левые части выражений (25) и (26). В резуль- 
тате получим: 


2 2 1 
(21) + (0) = т + а Хтиту). (28) 
Перемножим между собой левые и правые части. 
2 2 2 ыы 
(6%>) (01) ео =) (то лто >) . (29) 


Полученные выражения позволяют точно вычислить значения обобщенных масс т, и 
ту › . Для наглядной геометрической интерпретации совместного решения выражения (28) и (29) 
удобнее переписать в виде: 


— д: (7) (7) — 
оо = 4} СИ + С, 5, = В, (30) 


2 2 2 2 
- ао — ео: в ве — ео) (О, >) = (1) ] 
и 2 2 7 и 2 2 ь 
(0.2) (© 1) (0..2) (О) 


Две пары попарно симметричных корня системы получаются на основе определения точек 
пересечения гиперболы (первое уравнение) с прямой (второе уравнение) (рис.5). Оценим обра- 








где 


зующуюся погрешность при определении обобщенных масс, предположив, что с, =0. Для это- 
го рассмотрим характеристический полином системы после замены р= /®, т.е. для стационар- 
ных периодических движений будем считать, что обобщенные массы нам заданы. Рассмотрим го- 
дограф Л, (©) рик (о): ДА, (@) = (с; то (с ‚ -ту 26”) - (< ®,)”. Наблюдаемыми являют- 


ся резонансные частоты системы (рис.5). 
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а) 
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Рис.5. Геометрическая интерпретация вычисления обобщенных масс (а) 


2 
и смещения резонансных частот в зависимости от а (6) 


Как видно, по мере увеличения (с®),)? нижняя резонансная частота смещается в низко- 


0 
частотную область (от ет к 0. ), а верхняя - в высокочастотную (от Фе к 2. › ). Эти часто- 


ь 0 
ты являются близкими лишь при (с. )? > 0. Заметим, что в теории колебаний частоты Фь и 


0 
Фу называются парциальными, а частоты же 6 и 0.2 — собственными частотами системы. 


Что касается смещения резонансных частот, вызванных вариациями матрицы демпфирования, то 
эти смещения малы. 

Изложенные математические модели и методика идентификации параметров матриц 
инерционных коэффициентов и матриц демпфирования позволяют проанализировать свойства 
динамических моделей подсистем металлорежущих станков при варьировании точки равновесия 
системы. 

Особенности динамических моделей подсистемы заготовки. Все приведенные данные и 
выводы относятся к подсистеме режущего инструмента. Исследования показывают, что они спра- 
ведливы и для подсистемы обрабатываемой заготовки. При рассмотрении подсистемы обрабаты- 
ваемой заготовки в некоторой фиксированной точке контакта заготовки с инструментом имеет 


место не только равенство 77%; — 77% › , но и дополнительно о —- 2 И Со, —С02. В этом 
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случае в пространстве У любая ортогональная система координат является главной. Однако при 
рассмотрении подсистемы обрабатываемой заготовки имеется дополнительная сложность, свя- 


занная с тем, что жесткость Со; =Су› является функцией координаты перемещения суппорта 


вдоль оси вращения заготовки. Достаточно проанализировать распределение обобщенных масс и 
жесткостей, обусловленных силой Е, = 200,0 кг, приведенное на рис.6, для различных точек кон- 


такта инструмента с заготовкой. 





со 10°, кг/мм 


























ту, КГ.С2/мм 

10,5 

9,0 0,15 

У: 

6,0 0,10 

4,5 

3,0 0,05 

1,5 

0,0 


150,0 300,0 450,0 600,0 750,0 мм 


Рис.6. Изменения приведенных коэффициентов жесткости и массы в зависимости 


от координаты 1 вдоль оси вращения заготовки 


Анализ идентифицированных параметров инерционных коэффициентов и коэффици- 
ентов диссипации. Проанализируем основные закономерности изменения матриц т, й в зави- 
симости от смещения точки равновесия. Вначале проанализируем матрицы 7 на примере под- 
системы режущего инструмента токарно-револьверного станка 1К62. На рис.7 приведены измене- 
ния инерционных коэффициентов в зависимости от модуля смещения точки равновесия 


Уо =./ (1)? + (у)? . Там же даны значения дисперсии, рассматриваемой в относительных едини- 
цах к математическому ожиданию. Необходимо обратить внимание на два обстоятельства. Во- 
первых, значения обобщенных масс по направлениям ^, и у практически равны между собой. 


Причем их изменения по мере смещения точки равновесия системы заметны лишь в окрестности 
нулевого смещения, т.е. при отсутствии постоянной составляющей внешней силы. При этом соб- 
ственные частоты колебаний в вариациях относительно точки равновесия заметно меняются. Во- 
вторых, дисперсия при оценивании обобщенных масс уменьшается по мере смещения точки рав- 
новесия от нулевого значения, т.е. по мере увеличения внешней силы, действующей на вершину 
режущего инструмента. Практическое равенство обобщенных масс свидетельствует о том, что 
матрицу масс в пространстве Х можно считать диагональной. Это вытекает непосредственно из 
анализа уравнения (4). 
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0,8 














0,15 


0,10 


0,05 





0,0 








Рис.7. Изменение обобщенных масс То и Ту > (а) и дисперсии 


при оценивании обобщенных масс (6) при смещении точки равновесия системы 


Демпфирование является избирательным по направлениям колебательных смещений. 
Увеличение жесткости приводит к уменьшению демпфирования. Именно поэтому демпфирование, 
как и упругость, можно моделировать на основе рассмотрения эллипса демпфирования, т.е. ана- 
логично рассмотрению жесткости деформационных смещений. На рис.8 приведены характеристи- 


ки изменения коэффициентов демпфирования по и йо2 по направлениям у, и у, в зависимо- 


сти от 7, =\/ (у, + (у,)* . 


Йо, Кг. С/мм 


0,8 





Рис.8. Изменение обобщенных коэффициентов демпфирования йо и йо2 


смещения модуля точки равновесия системы 
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Заключение. Обобщая приведенный материал, можно заключить, что для подсистемы режущего 
инструмента в плоскости, нормальной к идеальной оси вращения заготовки в рассматриваемом 
частотном диапазоне, справедлива схематизация динамической системы без связи, формируемой 
процессом резания (см. рис.1). Важно подчеркнуть, что матрицы Й и с связаны между собой ли- 
нейным преобразованием, определяемым матрицами угловых коэффициентов. В этом случае име- 
ется возможность одновременного приведения к каноническому виду трех квадратичных форм: 
кинетической, потенциальной энергии и диссипативной функции Релея, т.е. существуют единые 
нормальные координаты, в нашем случае, — это координаты пространства Л. Указанное свойст- 
во, справедливое для рассматриваемых подсистем лишь при моделировании первых форм коле- 
баний, основано на неизменности ориентации эллипсов жесткости и демпфирования. 

Основные результаты исследований получены при финансовой поддержке РФФИ по про- 
екту 07-09-90000. 
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